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Abstract. We consider a nonlinear model for electrical conduction in biological 
tissues. The nonlinearity appears in the interface condition prescribed on the cell 
membrane. 

The purpose of this paper is proving asymptotic convergence for large times to 
a periodic solution when time-periodic boundary data are assigned. The novelty 
here is that we allow the nonlinearity to be noncoercive. We consider both the 
homogenized and the non-homogenized version of the problem. 
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1. Introduction 

We study here a problem arising in electrical conduction in biological tissues with the 
purpose of obtaining some useful results for applications in electrical tomography, see 
0|, [5], 0, 0, 0: 0, m, HU, [E], [13!: [H]- Our interest in this framework is mo¬ 
tivated by the fact that composite materials have widespread applications in science 
and technology and, for this reason, they have been extensively studied especially 
using homogenization techniques. 

From a physical point of view our problem consists in the study of the electrical 
currents crossing a living tissue when an electrical potential is applied at the boundary 
(see cm, [19] , [22], [25], [29j). Here the living tissue is regarded as a composite 
periodic domain made of extracellular and intracellular materials (both assumed to be 
conductive, possibly with different conductivities) separated by a lipidic membrane 
which experiments prove to exhibit both conductive (due to ionic channels in the 
membrane) and capacitive behavior. The periodic microstructure calls for the use of 
an homogenization technique. Among the wide literature on this topic, we recall for 
instance P, [2], [3], [IS], [E], [IS], [2D], [21], [28], [30], pH], [32], [33], [3]. Asa result 
of the homogenization procedure we obtain a system of partial differential equations 
satisfied by the macroscopic electrical potential u , which is the limit of the electrical 
potential u e in the tissue as e (the characteristic length of the cell) tends to zero. 
Different scalings may appear in this homogenization procedure and they are studied 
in |9j and [13]. We study here further developments of the model proposed in [4], 
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0 , 0 , 0 , n, pn, [13] . where the magnetic held is neglected (as suggested by 
experimental evidence) and the potential u £ is assumed to satisfy an elliptic equation 
both in the intracellular and in the extracellular domain (see, (12.111 below) while, on 
the membranes it satisfies the equation 



where [u £ ] denotes the jump of the potential across the membranes and cr £ X7u £ • v £ 
is the current crossing the membranes. From a mathematical point of view a big 
difference does exist between the case of linear / and the nonlinear case, as already 
pointed out in ra and P3- 

At least in the linear case, the asymptotic behavior of the potentials u £ and u is 
crucial in order to validate the phenomenological model employed in bioimpedance 
tomography devices, which currently relies on the use of complex elliptic equations, 
see P2]-P2]. 

Motivated by the previous considerations, in p3] and in this paper we investigate the 
behavior as t —> +oo of the nonlinear problem introduced in [13] . 

In 114) . we proved that, if periodic boundary data are assigned and / is coercive in 
the following sense 


feC\n), f'(s)>K> 0, \/seR, 


( 1 . 1 ) 


for a suitable k > 0, then the solution of the e-problem converges as £ —>■ +oo 
to a periodic function solving a suitable system of equations. In that case such a 
convergence was proved to be exponential. A similar asymptotic exponential behavior 
was proved for the solution of the homogenized problem. Similar results in different 
frameworks can be found in [23], m, ra, EH- 

It is important to note that in [10j-[12] our approach was based on eigenvalue esti¬ 
mates which made it possible to keep into account (as far as the asymptotic rate of 
convergence is concerned) both the dissipative properties of the intra/extra cellular 
phases and the dissipative properties of the membranes. 

Instead, in the nonlinear but coercive case, we proceed by exploiting the coercivity 
of /, hence the electrical properties of the intra/extra cellular phases do not appear 
in the rate of convergence. 

If / is not coercive it must be assumed to be monotone increasing and we proceed 
via a Liapunov-style technique so that the rate of convergence is not quantified. 

The paper is organized as follows: in Section [2] we present the geometrical setting 
and the nonlinear differential model governing our problem at the microscale e. In 
Section [3] we prove the decay in time of the solution of the microscopic problem. 
Finally, in Section H] we prove the decay in time of the solution of the macroscopic 
(or homogenized) problem, providing also the differential system satisfied by such 
asymptotic limit. 




ASYMPTOTIC ... 


3 


2. Preliminaries 

Let ft be an open bounded subset of R N . In the sequel 7 or 7 will denote constants 
which may vary from line to line and which depend on the characteristic parameters 
of the problem, but which are independent of the quantities tending to zero, such as 
£, 5 and so on, unless explicitly specified. 


2.1. The geometrical setting. The typical geometry we have in mind is depicted 
in Figure [TJ In order to be more specific, assume N > 2 and let us introduce a 
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FIGURE 1. On the left: an example of admissible periodic unit cell Y = 

E\ UE 2 UT in R 2 . Here E\ is the shaded region and E is its boundary. The 
remaining part of Y (the white region) is E 2 . On the right: the corresponding 
domain ft = ft\ U i2f U E e . Here I7f is the shaded region and T £ is its 
boundary. The remaining part of ft (the white region) is I7|- 

periodic open subset E of R N , so that E + z = E for all ;j G Z N . For all e > 0 define 
i?} = ft fl eE, i ?2 = ft \ eE. We assume that ft, E have regular boundary, say of 
class C°° for the sake of simplicity, and dist(r £ , dft) > 'ye, where r £ = dft{. We also 
employ the notation Y = (0,1)^, and Ei — E C\Y , E 2 — Y \ E, T = dE 0 Y. As a 
simplifying assumption, we stipulate that E\ is a connected smooth subset of Y such 
that dist(r, dY) > 0. We denote by v the normal unit vector to E pointing into E 2 , 
so that is e (x) = z/(e _ 1 a:). 

For later use, we introduce also the conductivity 
J 0-1 if y G E u 

a(y) = < and a 0 = \Ex\ai + \E 2 \<t 2 , 

1 ^0-2 if y G E 2 , 

where o"i,o- 2 are positive constants, and we also set cr £ (x) = a^e^x). Moreover, let 
us set 

£^(y) {u : Y \ r —> R \ u\e x 6 C k (Ei) , u\e 2 € C k (E 2 ) , and u is Y — periodic} , 
for every 0 < k < + 00 , and 

X#(Y) := {u G L 2 (Y) I u\e x £ H l {E \), u\e 2 £ H l {E 2 ), and u is Y — periodic} . 
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More generally, the subscript # in the definition of a function space will denote 
periodicity with respect to the first domain, in such a way that the extended function 


remains (locally) in the same space. 

We set also 

X 1 {fi £ ) := { U £ L 2 (f}) | U\Q e G U\Qe G 2 )} • 

We note that, if u G Xl (Y) then the traces of u\E t on r, for i — 1,2, belong to 
H 1 ^ 2 (r), as well as u G A’ 1 (i7 £ ) implies that the traces of uiq? on r £ , for % — 1,2, 
belong to H l / 2 {r e ). 

2.2. Statement of the problem. We write down the model problem: 

— div(< 7 e Vu £ ) = 0 , in (i?f U 12|) x (0, T); (2.1) 

[a £ Vu £ ■ v e \ = 0 , on r £ x (0, T); (2.2) 

+ =ff£ V?i £ ^ £ , on r £ x (0, T); (2.3) 

[u e ](z, 0) = 5 £ (a:), on T £ ; (2.4) 

u e (x) = \l/(a;,t), on dft x (0, T), (2.5) 


where <r e is defined in the previous subsection and a > 0 is a constant. We note that, 
by the definition already given in the previous section, u e is the normal unit vector 
to r £ pointing into i?f. Since u £ is not in general continuous across r e we set 

:= trace of u £ \q^ on r £ x (0, T); m( 2) := trace of u £ \qs on r £ x (0 ,T). 

Indeed we refer conventionally to as to the interior domain, and to i?f as to the 
outer domain. We also denote 


Similar conventions are employed for other quantities, for example in (12.21) . In this 
framework we will assume that 

i) S £ G H 1 / 2 (r £ ), ii) J S 2 (x ) da < 7 £ , ( 2 . 6 ) 

where the second assumption in (12.61) is needed in order that the solution of system 


(12.11) (12.51) satisfies the classical energy inequality, (see (3.1) in [14]). 

Moreover, / : R —» R satisfies 

feC\R), (2.7) 

/ is a strictly monotone increasing function, ( 2 - 8 ) 

/(0)=0, (2.9) 

f'(s) > S 0 , for a suitable h 0 > 0 and V|s| sufficiently large. ( 2 - 10 ) 

The previous assumptions imply also 

f(s)s > Ai s 2 — A 2 |s| , for some constants Ai > 0 and A 2 > 0 . ( 2 - 11 ) 
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Notice that the results presented in this paper hold also in a more general case, namely 
if we replace condition (12.101) with the assumption that f^ 1 is uniformly continuous 
in R. for example when f(s) = s + sins. 

Finally, T : i? x R —> R is a function satisfying the following assumptions 

i) <P6LL(«;» 2 («)); 

«) ( 2 . 12 ) 

Hi) 4/(a;, •) is 1-periodic for a.e. x G i?. 

Existence and uniqueness for problem (12.1[) (12.51) has been proved in [8]. Moreover, by 
m Lemma 4.1 and Remark 4.2] it follows that the solution u £ G C°((0, T]; fh 1 (l7 e )) 
and [u £ ] G C°((0, T]; L 2 (F e )), uniformly with respect to e, and [u £ ] G C°([0, T]; L 2 (r ,£ )), 
but with non uniform estimates. 

3. Asymptotic convergence to a periodic solution of the £-problem 

The purpose of this section is to prove the asymptotic convergence of the solution of 
problem (12.ID — f 1 2.5 [) to a periodic function uf when t —> +oo. The function uf is, in 


turn, a solution of the system 

— div(<7 e Vuf) = 0 , in (i?^ U i?|) x R ; (3.1) 

[cr £ X7uf ■ u £ ] — 0, on r £ x R; (3.2) 

~[«f] + / = (ff'Vuf ■ V .) , on r' X «; (3.3) 

uf(x,t) = 'i/(x,t), on dfl x J2; (3.4) 

uf(x, •) is 1-periodic, ini?. (3.5) 


Indeed, this problem is derived from (12.1j) — (12.51) replacing equation (12.41) with (13.5j) . 
The rigorous dehnition of weak solution of (13. 1 j) — (13.5D is standard (see for instance 
P3, Dehnition 4.13] 

As a first step we will prove the following result. 

Proposition 3.1. Under the assumptions (12.71) - (12.101) and (12.121) . problem (13.11) — 
(I3.5D admits a solution uf G C^([0,1]; A 1 (i? e )). 

Proof. For 5 > 0, let us denote by fs(s) := /(s) + Ss, for every s G R, and consider 


the problem 

— div(cr £ Vu^ 5 ) = 0 , in (i?^ U i?|) x R; (3.6) 

[o~ £ Vuf s ■ uf\ — 0 , on r £ x R ; (3.7) 

+ fs (^4^ j = a£Vu ts • > on r x R > ( 3 - 8 ) 

uf s (x, t) = T(a;, t) , on <9i? x R; (3.9) 

uf s (x, •) is 1-periodic, ini?. (3.10) 
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For any positive e and S, the previous problem admits a unique time-periodic solution 
because of the results already proved in [14] . 

On the other hand, multiplying equation (13.61) by uf s — \&, integrating by parts on 
i? x [0,1], using the periodicity and taking into account equations (I3.7I) - (I3.9D . we get 

// Y^ u * s ^ dxdt + J J ( ~T ~) ^ dadt - J j yl v *l 2 dxdt. (3.11) 

on or*'' on 

Finally, using (12.111) we obtain 

l l l 

If \\Vu* 5 \ 2 dxdt+ J j ^[ufs\ 2 dadt < J j \Vd>\ 2 dx dt+^-X 2 2 \r £ \ . (3.12) 

on on on 

Multiplying now equation (13.61) by uf St — d> t , integrating by parts on i? x [0,1], using 
the periodicity and taking into account equations (I3.7D - 03.9D . we get 



o n on on 

i i 


-// yl VvI/ | 2dxdt +^ A 2l r£ l + /f y|V^| 2 dxdO (3.13) 
on on 

where we used (j3 . 12[) . Notice that the second integral on the left-hand side is equal 
to zero by periodicity and trivial integration. Hence 

l l l 

Y e f J l u t,5,t} 2 dt < J J y|V'h t | 2 da:dt + J J dx dt + ^-\ 2 2 \r e \. 

or® o n on 

(3.14) 

Inequalities (13.12j) and (13.141) , for e > 0 hxed, yield the weak convergence of uf g and 
Vuf s in L 2 {Qf x (0,1)), % — 1,2, and respectively the strong convergence of [uf s ] 
in L 2 (r e x (0,1)), for <5 —>■ 0. Since all the functions uf s are 1-periodic, denoting as 
usual with uf the limit of uf s we have that the same periodicity holds true for uf. 
Moreover we can pass to the limit, as 8 —» 0, in the weak formulation of problem 
(I3.6D (13. 10D . thus obtaining that uf is a 1-periodic solution of problem (13. ID — (I3.5D . 
under the assumptions (I2.7D - 02.1OD and (12.121) . 
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Differentiating formally with respect to t (l3.6IMI3.9p . multiplying the first equation 
thus obtained by (uf 5t — \Eq) and finally integrating by parts, we obtain 

i 1 

/ f Y\ Vu ts/ dxdt ^ f J y|V* t | 2 <Lrd*. 

on on 

Since the estimates above are uniform in 5, we have that uf belongs to the class 
claimed in the statement. □ 


Given £ > 0, it remains to prove the asymptotic convergence of the solution u £ of 
(I2.ip ~ 02.5p to uf, for t —> +oo. 

Theorem 3.2. Let £ > 0 be fixed and let u £ be the solution of problem (12.11) - 02.5p . 
Then, for t —> +oo, u £ —> uf in the following sense: 


lim | \u £ (-,t) -uf(-,t)\\ L 2 {{2) = 0; 

t — H-oo 

lim || Vu £ (-, t) - \7uf(-,t)\\ L 2( n) = 0; 

£—>•+00 

H[“'](' ,t ) - tiller*) =0. 


(3.15) 

(3.16) 

(3.17) 


Proof. Setting r £ := uf — u £ , we obtain that r £ satisfies 

— div(cr £ Vr e ) = 0 , in (i+ U i?^) x (0, +oo); (3.18) 

[<7 £ Vr £ • v £ ] — 0 , on T £ x (0, +oo); (3.19) 

qj n \v \ 

--QfVe] + 9e(x, ty-j- = cr £ Vr e • u £ , on T £ x (0, +oo); (3.20) 

Ve}{x, 0) = [uf(x, 0)] — S E (x) =: S £ (x), on D; (3.21) 

on <9i? x (0, +oo); (3.22) 

if [u £ ](x,t) = [uf}(x,t), 

if [u £ ](x,t) f [uf](x,t), 

so that g £ (x,t ) > 0, and S e (x) still satishes assumption ii) in (12.61) because of (13.121) 
and (13.141) . Multiplying equation (13.181) by r £ and integrating by parts we have 

j cr e \Vr £ \ 2 dx+ ^ J [r £ fi[r £ ]da + J 9c ^ ^ [r £ ] 2 da = 0. (3.23) 

n r e 

Equation (13.231) implies that the function (4^ J re [r £ (x, t)} 2 da is a positive, decreas¬ 
ing function of t; hence, it tends to a limit value r £ > 0 as t —> +oo. We claim that 


r £ (x) = 0, 


where 


u 


f — Om) 


9e(x,t) := < 




[ {x,t ) - ^(z,t) 
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the value r £ must be zero. Otherwise, for every t > 0, fire [r e (x,t)} 2 da > r £ > 0. 
Ou the other hand, setting (t) := {i 6 P : [ r £ (x,t )] 2 < 0 ^p e \ }, we have that 


y- J [r £ (x,t)] 2 da > y , Vt > 0 . 

r e \ri{t) 


(3.24) 


Indeed, by definition, 

ol (, , Sl9 , a 
r £ < ~ [r e (x,t)] 2 da = - 
el e 


a 


[r £ (x,t)} 2 da + - / [r e (x, t)] 2 da 


r e 


< 


a 


r e \rf e (t) 

r / ., 9 a rye . „ . a 




r s \rs(t) 


e 2a\r 


[ry(x,t)] 2 da + , 


r s \rf(t) 


which implies (13.2411 . Moreover, we have that, on r £ \ (t), g £ (x,t) > x > 0, where 

X is a suitable positive constant depending only on (r e ,£,at,\r e \) (this last result 
follows from assumption (l2.8fM2.10p ). Hence, using (I3.23|) . it follows 


cl I a 
dt 1 2e 


[r e (x,t)} 2 da < - 


9e(x,t) 


r £ (x,t)] 2 da 


r e 




< -x 


~[ r £ ( x , t)} 2 da < —yy-y < 0. (3.25) 

£ ZOL 


r e \ri(t) 


Inequality (13.2511 clearly contradicts the asymptotic convergence in t of the function 
t 1 y ~ f re [r £ (x, t)] 2 da, hence 


a 


lim — / [r £ (x, t)] 2 da = 0 . 

£—>■+00 £ J 


(3.26) 


r e 


In particular, this gives (13.171) . Integrating (j3.23jl in [t, oo) and taking into account 
(13.261) . we get 

+oo 



a 


cx e |Vr e | 2 d:rdt < — / [r e (x, t)] 2 da , 


t ^ 


re 


which implies 


+oo 

lim / / a £ \\7r e \ 2 dx dt = 0 

t D 


(3.27) 


(3.28) 


Condition (13.28jl guarantees that for every positive rj there exists a t(rj) > 0, such 
that 

+oo 

a e |Vr e | 2 dxdt < rj, 



f(v) 


n 
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which, in turn implies that, for every natural number n, there exists a t n G ( t{j) ) + 
n, t(rj) + (n + 1)), such that 


a £ \Vr e (x,t n )\ 2 dx < 17 . 


(3.29) 


a 


Now, we multiply (13.181) by r £ji and integrate in i2, so that 


n 


a £ X7r £ Vr £tt {x, t) dx + ^ ^[r £jt (x, t )] 2 da + J 

pe pe 


r £ (x,t)} [r E:t (x,t)}da = 0, 


which implies 


a £ Vry'Vry t dx < 


9l(x,t) r l2 


2ae 


[r £ ] da. 


(3.30) 

(3.31) 


n r e 

Moreover, integrating (I3.3ip in [t n ,t*] with t* G [t n ,t n + 2] and using (13.29[) . we have 


sup 

tG [tn ,?n +2] 


a c 

T 


?? L 2 

|Vr £ (x,t)| 2 dx | < ^ sup 


a 


Vi? 


2 t€[tn,+Oo) l ^ 


[r £ (x, t )} 2 da , \/n e N . 


r e 


Since t n+1 — t n < 2, the intervals of the form [t n ,t n + 2], when n varies in N , are 
overlapping; hence, we obtain 


sup 

££ [t-hljH-oo) 


a 


17 L 


a 


— |Vr e (a;,t)| dx < - + — sup - [r e (x,t)] da . (3.32) 


\n 


2 a te[t,+oo) \ £ 


r e 


Because of (13.261) the integral in the right-hand side of (13. 3 2D can be made smaller 


-1 


than 2 \ ■ provided t is chosen sufficiently large in dependence of ?/. This means 

that 


sup 

££ [£-l-l,-l-oo) 


a' 


— \S7r £ (x,t)\ 2 dx | <77, 




so that 


lim / cr £ \X7r £ (x, t )\ 2 dx = 0 

t —^-|-oo J 
Q 


(3.33) 


(3.34) 


I 11 particular, this gives (I3.16D . Finally, Poincare’s inequality together with (13.261) 
and (I3.34D yield 


lim / \r £ (x, t )\ 2 dx = 0 , 

t —(“OO J 

f? 


which gives (I3.15D . 


(3.35) 

□ 


Remark 3.3. More in general, the previous procedure allows us to prove that solutions 
of (I2.1D - (I2.5D having different initial data satisfying (12.6D but the same boundary 
condition tend asymptotically one to the other (such convergence being exponential 
if / is coercive in the sense of (II.ID ). □ 
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Remark 3.4. Observe that, thanks to previous remark, Theorem 13.21 implies unique¬ 
ness of the periodic solution of problem (I3.1IHI3.5D in C^([0,1]; fh 1 (i7 e )). □ 


4. Asymptotic decay of the solution of the homogenized problem 


The aim of this section is to prove asymptotic decay of the solution of the homogenized 
problem. To this purpose, let (u, u 1 ) G L 2 ((0, T); // 1 (i7)) x L 2 (i7 x (0, T); X#(Y)) be 
the two-scale limit of the solution u e of problem (I2.1IMI2.5I1 . where the initial data S £ 
satisfies the additional condition that S £ /e two-scale converges in L 2 (i?; L 2 (T)) to a 
function Si such that Si(x, •) = S\r(x, •) for some S G C(i7; C^(T')), and 

hr/ (§) w<n= 

r e nr 



S i(x, y ) da; dcr. 


(4.1) 


We recall that, under these assumptions, by m Theorem 2.1], the pair ( u,u x ) is the 
weak solution of the two-scale problem 


a 0 Vu + Jcr'VyU 1 dy\ 

= 0, 

in i? x (0, T); 

(4.2) 

Y / 

— diVy(aVw + aVyU 1 ) 

= 0, 

in i? x (Ai U E 2 ) x (0, T); 

(4.3) 

[ff(V« + VyW 1 ) ' P] 

= 0, 

on 17 x h x (0, T); 

(4.4) 

a ^i ul ] + f {[u 1 ]) 

= a(Vu + V.ytt 1 ) 

• p , on i? x r x (0, T); 

(4.5) 

[m 1 ] (a;, 2 /, 0) 

= Si(x,y) , 

on 17 x h; 

(4.6) 

u(x,t ) 

= V(x,t ), 

on <9i? x (0, T); 

(4.7) 


in the sense of the following definition. 


Definition 4.1. A pair (u, u 1 ) g7 2 (( 0, T); i/ 1 (i7)) x L 2 (i7 x (0, T); is a weak 

solution of (14.2(1(14.7(1 if 



o nr nr 


for any function (j) G C°(0, T; //,)(]?)) and any function <f> G C°([0, T\; L 2 (i 7; X#(Y))) 
with [<f> t ] G C°([0, T]; L 2 (i? x T)) which vanishes at t = T. 

Moreover, u satisfies the boundary condition on <9i? x [0, T] in the trace sense (i.e. 
u(x, t) = v k(x, t ) a.e. on <9i? x (0, T )) and u 1 is periodic in Y and has zero mean value 
in Y for a.e. (x, t) G i? x (0,T) (see p3] Definition 5.1]). 
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For later use, let us define 


M 1 ^)) + H^ 1 ||c°([0,l];L 2 (17xY)) 

+ II Vyh ||c°([0,l];L 2 (i?xY)) + II [h ] ||e°([0,l];L 2 (12xr)) , (4-9) 

where (/r,/? 1 ) G C°([0, T]\H l (R)) x C°([0, T]; L 2 (J7; X#(Y))), and 

\\\\(h, ^ 1 )|||| := ||^||^ 1 (J7) + ||/l 1 ||L 2 (ftxY) + II V y h 1 || jfy 2( fi xY) + II [/l 1 ]||L 2 (12xr) , (4-10) 

where (hji 1 ) G i/ 1 (i7) x L 2 (i?; X#(Y)). 

As in the previous section, first we prove that there exists a time-periodic weak 
solution of the two-scale problem 


a 0 Vu # + J 

aVyU 1 ’# dy J 

= 0, 

in 

R x R\ 


(4,11) 

Y / 

— div y (a\7u* + a\7 y u l, rt) 

X 

.S 

o 

(E 

i U E 2 ) x 

R 

(4.12) 

[a(Vw # 

+ V,v>*) ■ „] 

= 0, 

on 

X 

Lh 

X 

R 

(4.13) 

d r i, 

“5 [ “ 

*} + / flu 1 "* 1 ]) 

= cr(V'U # + VyU 1 ’#) • ZZ 

,on 

X 

X 

R 

(4.14) 


\v}*}(x,y,-) 

is 1-periodic, 

on 

Rxr- 


(4.15) 


U*(x, t) 

= d>(x,t ), 

on 

OR x R; 


(4.16) 


in the sense of the following definition. 

Definition 4.2. A pair (v # , v 1 ’#) G C^QO, 1]; Ff 1 (i?)) x C#([0,1]; L 2 (i?; X#(Y)) with 
[Vf’^| G L#(0, 1; L 2 (i? x J 1 )) is a time-periodic weak solution (with period 1) of 

(imu- mri]) if 


J J I cr (Vv # (x,f) + \' y y 1 ’*(x,y,t)) (V<j>(x,t)+ V y ®(x,y,t)) dxdydt 

R R Y 

+ J J J f([v 1,# {x,y,t)])[$(x,y,t)]dx dadf 

r n r 

~ a J f f[ yl ’*( X,y,t ^^®( X,y,t ^ <h:d<Tdt = 0 ( 4 ' 17 ) 

r a r 

for every (0, <f>) E C° c (R; H^R))xC° c (R-, L 2 (R; X^Y)), [$ t ] G L 2 (R; L\Q x T)) and 
v 1,7 ^ has zero mean value in Y for a.e. (x, t) G R x R and v# satisfies (14.1611 in the 
trace sense (see HI Definition 5.7]). 

Remark 4.3. We note that by a standard approximation of periodic testing functions 
with functions compactly supported in a period, the weak formulation (14.1711 can be 
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equivalently rewritten as 


J J J a (Vv#(x,t)+ X7 y v 1 ' # (x,y,t)) (V<f>(x,t)+ V y $(x,y,t)) dxdydt 

0 17 Y 

1 

+ /// f([v 1,# (x,y,t)])[®(x,y,t)]dx da dt 
oar 

i 

-a// J[v 1 '*(x,y,t)\^[<S>(x,y,t)]dxdadt = 0 
oar 

for every (0, $) £ C° # {[ 0, x C%([ 0,1]; L\Q- X#(Y)), [$ t ] £ L%( 0,1; L\Q x 

r)Y Hence, when it is more convenient, we replace compactly supported testing 
functions with 1-periodic testing functions. □ 

Proposition 4.4. Under the assumptions (12.711 - (12.1011 and (12.121) . problem (14.1111 - 
(14.1611 admits a 1-periodic in time solution. 

Proof. For 5 > 0, let us denote by fs(s) := f(s ) + Ss, for every s £ R, and consider 
the problem 


- div I a 0 Vuf + / aV v uf v dy I = 0 , 


6 #. 


in 17 x R] 


Y 


(4.18) 


dw y (aVuf + aVyu]’*) = 0 , 
{a{X7uf + V y u]’*}u] = 0, 


in i? x (Ei U E 2 ) x R ; 


(4.19) 


on 17 x r x R; 


(4.20) 


a ^[ u l’*} + fs ([«J’ # ]) = <x(Wuf + Vyu]’*) -v, on R x T x R] 

(4.21) 

on 17 x F ; (4.22) 


ul’*](x, y, ■) is 1-periodic, 


uf(x,t) = d>(x,t), 


on <9l7 x R ; 


(4.23) 


where has zero mean value on Y for a.e. (x, t) £ 17 x R. 

Since fs has a strictly positive derivative on R , by the results proved in m Section 
5], a unique periodic solution (uf of problem (I4.18D - (14.231) does exist, i.e. 
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uf, u l'#) satisfies 



a (Vuf(x,t) + X7 y u L s ,# (x,y,t )) (V0(x,t) + V y ®(x,y,t)) dxdydt 


A,#t 


0 Q Y 



fsi[u\*{x , y , t)])[$(x, y, t)} dx da d t 


o si r 


a 



u]’*(x,y,t)]^[<f>(x,y,t)]dxdadt = 0, (4.24) 


d 


o a r 


for every (0, $) G C%{[ 0,1]; Hl(Q)) x C%([ 0,1]; L\Q- Aj(r)), [4> t ] G L\ (0,1; L 2 (f2 x 
Z 1 )) (recall Remark 14.31) . Moreover uf# has zero mean value in Y for a.e. (x,t) G 
Q x R and uf satisfies (14.231) in the trace sense. By (14.241) we get that (uf,uf*) 
satisfies an energy estimate, easily obtained replacing (0, 4>) with (uf — 40 w^), 
which implies 



Vuf + V y uf*\ 2 dxdydt + 



fs([ui*])[u]’ # ] dx da d t 


o n y 


oar 



a 


V\R | 2 dxdydt , (4.25) 


o a y 


where we take into account 

i 

i# 


9 r.dd 


KflK*]® = ^ / ^K # ] 2 d« = o, 


(4.26) 


0 


1 44 

which is a consequence of the periodicity of u$. 

From (14.251) . working as done in (13.Ill) - (13. 121) of Section [3] and taking into account 
( 12 . 111 ) we get 



a\Vuf + \/yuf^\ 2 dx dy dt + 
o a y oar 

where 7 is a constant depending on Ai, A 2 , \T\ and the H 1 -norm of 4C 

Replacing (0, 4>) in (I4.24p with (uf t — 4by (14.271) . (12.111) and taking into 

account the fact that 


1,#|2 



Ai[ m 5 ,# ] 2 da dt < 7 , (4.27) 


(Vuf + V 9 «7)(ViiJ t + V„ujf) dt = I / 0 1 V “f + Vdl'T dt = 0 


. 1 1 ~tj~ \ / V7 ,1 A 


1 ,# 


d 


# 




0 


0 
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and, denoting by F$ a primitive of f$, 

1 

J J J fs{[u 5 ,# })\u 5 f]dxdadt = 

o a r 

because of the periodicity, we get 

i 



aft (Id’*]) 

dt 


dx da dt = 0 , 


o a r 


a / [u]f} 2 dxdadt < 7 , 


(4.28) 


oar 


where, again 7 depends on Ai, A 2 , \T\ and the H 1 -norms of 4' and 4+ From (|4.27|h 
we obtain 



0 a 
1 



uf 2 dxdt < 7 , 

(4.29) 

2 dx dy dt < 7 . 

(4.30) 


0 a y 


Indeed, 



y,t)\' 2 dy dx dt + J j |V+f | 2 dxdt 
0 a 


0 a y 

1 


<7-2 


,#| 2 , 



X7yu]’#(x, y, t) V+f ( x , t) dy dx dt 


#/ 


0 a y 
1 



=7-2 / / Vuf I / V y ^’ ff (x,i/,t)di/ I dx dt 


1,#/ 


0 a 
1 


\Y 



<7 + 2 / / | Vwf | 


|[^’ # (x,?/,t)]| dcr I dx dt 


0 a 

1 1 

#12 



<7 + - / / |Vnf | 2 dxdt + 8 |-T| [-+5 (x, y, t )] 2 da dx dt 



0 12 
1 


0 D r 



| Vwf(x,t )| 2 dxdt + 7 . 


0 a 


(4.31) 
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In order to be able to pass to the limit 5 —> 0 we need a formulation with vanishing 

JX _LL J± 

boundary data. To this purpose we set vf = uj — T; clearly vj satisfies 


TjVuf + J aW y u]'*dy \ =div(a 0 VT) , 

in 17 x R 

Y J 

(4.32) 

- div, ,(aVvf + aV y u]'*) = 0 , 

in 1? x (Ei U E 2 ) x R 
(4.33) 

[a(Vuf + V y u)'*}u\ = — [aVT • v\ , 

on 1? x r x i? 
(4.34) 

+ fs (K’ # ]) = ^ 2 (Vuf + V,/uJ’ # ) 

• v + a 2 VT • i/, on 12 x f x ii 

(4.35) 

[u^’ # ](x, i/, •) is 1-periodic, 

on 17 x T 


(4.36) 

vf(x, t) = 0 , on dft x R , 

(4.37) 


or, in the weak form, 



a (Vuf (.x , t) + V y u^(x, y, t )) (V0(x, t) + V y <h(x, y, t)) dx d y d t 


,#/ 


o a y 



f s ([u]' # (x, y, t)])[$(x, y, t)} dx da d t 


o si r 
i 



a I I I [u)'*(x,y,t)\^\§(x,y,t)} dxdadt 
o or 


d 




a o VT(x,t)V0(x,t) dxdt + / / [aVd>(x,t) ■ u]^ 1 ' ) (x,y,t) dx da dt 


o i? 


o i? r 



o‘ 2 V'h(x, f) • z/[<h(x, y, f)] dx da dt , (4.38) 


o nr 


for (0, <E>) as in Remark 14.31 At this point, (j4.27jl fj4.30p allow us to pass to the limit 
with respect to S in the weak formulation f|4.38p . thus proving that there exists a peri¬ 
odic (in time) pair of functions (u # , u 1 '#) G L#( 0,1; Hi (17)) xh| (0,1; L 2 (i7; X#(Y))) 
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such that u 1 ^ has zero mean value on Y, for a.e. (x,t) E i? x R, and ( v^,u 1,# ) sat¬ 
isfies the homogenized problem 



cr 



(Vu # + VyU 1 ’#) ■ X7(j)dxdydt + / / a(Vv* + V y w 1,# ) ■ V y &dxdydt 


0 17xY 


0 OxY 

i 



//[$] dx dcr dt — a 



[u 1 ’*]^] dx da dt 


o nr 
i 


o nr 

i 




cr o V\l/V0 dx dt + / / / [crV'P • dcr dt 


o n 


o nr 


+ 



cr 2 V^ • zz[$] dx da dt (4.39) 


o nr 


for every test function (0, $) G C^([0,1]; Hq (17)) xC^([0,1]; L 2 (i7; X#(Y)) with [<h 4 ] G 
L?L(0,1; L 2 (i? x T)), where we have taken into account that (12.71) and (I4.27p imply 


fs(bi'*]) = /([uj" # ])+5[%’ # ] /I, weakly in L 2 ((0, T) x Q x T £ ), when <5-^0. 

It remains to identify //. To this purpose, we follow the Minty monotone operators 
method. Let us consider a sequence of 1-periodic in time test functions ij)k(x,y,t) = 
<j>o(x,t) + 4>\ (x,y,t) + A0 2 (x,y,t), with 0£ E C°°(i1 x R), (p'l E C°°(i7 x R ; <£#(F)), 
02 G C] (i? x (0,1); £^(Y)), with 0 q(-, t) vanishing on <9i? for t G -R, 0q —» strongly 
in tf oc (R; H^(Q)) , ^ —>■ u ll# strongly in L 2 oc (R] L 2 (i7; A^(y))), [0j] —>■ |V’ # ] and 
[0* t ] —> [uj’ # ] strongly in L 2 ((0,1) x i? x T), i.e. 


l 

J f ~u 1 ,# (x,-,t)\\ 2 Hl{Ei) dtdx 

o i? 

i 

+ J J \\[<f>i(x,-,t)] - [u 1># (a;,-,i)]||| 2 (r) dtda; 

0 17 

1 

+ J J \\[&lt(xr,t)]-[u]’*(x,-,t)]\\ 2 L2{r) dtdx ^>- 0 , for k ->■ +oo, i = 1, 2. 

0 17 

Clearly, 0 q can be constructed by means of standard convolutions with regular kernels; 
instead, in order to construct 0^ we proceed as follows. Taking into account that, 
passing to the limit for 5 —> 0 in (I4.28p . we have 

[u]'*} E L 2 ((0,1) x n x r) , (4.40) 

by standard arguments we can approximate the jump [u l '^\ with a sequence of 1- 
periodic in time functions (j)\ E C°°(i 1 x T x R) such that 0^ —> [w 1 ’#] strongly in 
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L 2 (( 0,1) x i?; H l / 2 (r)) and 0^ — >• [i u\'*} strongly in L 2 ((0,1) x ft x I 1 ). Now, define 
( f>i as the (1-periodic in time) solution of the problem 


div„(«T(V^ + V„0S) =0, 

in (Ei U E 2 ) x i? x R; 

(4.41) 

[ff(V4 + V v *f) ■ v] = -[(TV* ■ v ], 

on r x f] x R; 

(4.42) 

W] = ft , 

on r x i? x R; 

(4.43) 


and is H-periodic with zero mean value on Y for (x,t) G i? x R. By 

Lemma 7.3 in m, it follows that <p\ G C°°(i 7 x [0,1]; £#(R)). Here, for the sake 
of simplicity, we work as if 4/ has enough regularity, otherwise we proceed with a 
standard regularization procedure also on T. Moreover, by [HI Lemma 5] applied to 
0i — u 1,# with P = div y (cr(V</>o — Vx # )) = 0 in E\ U E 2 , Q = [cx(V0o — Vn # )], and 
S = 4>i ~ [u 1,# ], we obtain 

l|01 — u1 ’ # \\l 2 ((0,1)xO-X^(Y)) < 7(l|01 — [ wl ’ # ]||L 2 ((O,l)xr2;ff 1 / 2 (r)) + l|V0o — Vu # ||L2((o,l)xfl)) • 

(4.44) 

Since the right-hand side of (14.44j) tends to zero for k —> +oo we obtain the desired 
approximation. 

Taking only into account the monotonicity assumption on /, the periodicity in time 
of (f )q and (j)\ and Remark 14.31 we calculate 


[ f <r(Vvf + X7 y ul’* - V0S - V y 0* - A V y 0 2 ) • (Vvf - V0S) dx d y d t 

0 QxY 
1 

+ [ f C r(Vnf + V y ^ # -V^-V^-AV s/ 0 2 )-(V^ # -V^-AV y 0 2 )dxd 1 /dt 

o nxY 

i 

+ « f f ^ ([%’ # ] - [01 + A0 2 ]) ([?4’ # ] - [01 + A0 2 ]) dx dcr dt 

o nxr 
i 

+ f f (fsdu]’*}) - fs ([01 + A0 2 ])) ([uj ,# ] - [01 + A0 2 ]) dx da dt 

o Qxr 

i 



a\Vvf + V 


y u 8 


h# 


- V0£ - V.^ - AV y 0 2 | 2 dx dy dt 


0 QxY 


+ f f (/«([«J’ # ]) - fs ([01 + A0 2 ])) ([uj’ # ] - [01 + A0 2 ]) dxdadt > 0 , (4.45) 

o Qxr 
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where we have taken into account that the time-periodicity of u 


i,# 

5 i 


0i and 0 2 implies 


a [ [ ^ ([^’ # ] - [01 + A0 2 ]) ([%’ # ] - [01 + A0 2 ]) dxdcrdt 

0 17 xf 


= f [ [ + A ^]) dxdadt 

o i?xr 

= | / (K’ # (x,!/, 1)] - [01 Or,?/, 1)]) dxdcr 

i7 xr 

/ (K’ # (x,i/,0)] - [^(x,i/,0)]) dx dcr 

17 xf 


0. 


Taking the function (xf" — 0 q — 0 1 — A0 2 ) as a test function ( 97 ,$) in (14.38)1 . 
inequality (I4.45jl can be rewritten as 


1 

[ [ cr(V<t>Q + V„0j + AV,0 2 ) • (Vnf - V0g) dx dy dt 

0 SlxY 
1 


- [ I a(V0^ + V^ + AV,0 2 )-(V^’ # -V y 0f-AV,0 2 )dxdydt 

0 i7xy 

1 

- a j f f J^[0t + A0 2 ] ([uj’ # ] - [0i + A0 2 ]) dx da dt 

0 17 l7xT 
1 


- f f fs([<t> 1 + A0 2 ]) ([«J’ # ] - [0? + A0 2 ]) dx da dt 

0 i7xr 


> J J ^oVT • (Wvf 

0 17 


V0q) dx dt 


/ / / [<w* • "M* 

0 17 r 


0i — A0 2 )dx dcr dt 


0 17 r 


cr 2 V\l/ • v(\u\’*] — [0i + A0 2 ]) dx dcr dt. (4.46) 
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Hence, passing to the limit as 6 —> 0 and using (14.281) . it follows 
a(V0* + Vy(f>\ + AV„0 2 ) • (Vu # - V0£) dxdydt 



0 fixY 
1 



ct(V0q + V y 0? + AV y 0 2 ) • (VX’ # - V„<# - AV„0 2 ) dxdydt 


0 12xY 


— a 



d 


— [0i + A0 2 ] ([w 1,# ] - [0t + A0 2 ]) dxdadt 


0 12 12xT 

1 



/ ([01 + A0 2 ]) ([w 1,# ] - [01 + A0 2 ]) dx dcr dt 


0 12x.T 

1 1 
> J J 0‘oV'h • (Vu # - V0g) dx dt - 

o a 



[aVT • zi](-u 1,# — 0 i — A 02) ( ' 1 ^ dxdadt 


o 12 r 



a 2 Vd> ■ z/([u 1,# ] — [0i + A0 2 ]) dx da dt. (4.47) 


oar 

Now, letting k —y +oo, we obtain 
i 

r r 

a(Vx # + VyU 1 ’* + AV y 0 2 ) • AV<,02 dx dy dt+ 



0 12 Y 
i 

a 



^ [n 1,# + A0 2 ] A [0 2 ] dx da dt + 


i 



/ ([w 1,# + A0 2 ])A[0 2 ] dx da dt 


o 12 r 


o 12 r 

i 


> 




[aVf • p]A 02 1 ' ) dx da dt + 
o 12 r oar 

Taking into account (I4.39j) with 0 = 0 and $ = 0 2 , (I4.48j) becomes 


a 2 VT • p[A 0 2 ] dx da dt. (4.48) 



A 2 / / / aV y 0 2 • Vy0 2 dx dy dt + aA 2 



dt 


[0 2 ] [02] dx da dt 


0 12 Y 

1 


o 12 r 


-A 




y[0 2 ] dx da dt + A / / / / ([u 1,# + A0 2 ]) [02] dxdadt > 0 . (4.49) 


o 12 r 


o 12 r 
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Assuming firstly that A > 0 and then A < 0, dividing by A the previous equation and 
then letting A —» 0, we obtain 

i i 



/i [ 02 ] da; da dt 


o i 2 r 


/ / /•^([' ul ’ # ])[^ 2 ] dxdadt, 
o i 2 r 


which gives 

h = / ([« 1,# ]) • (4-50) 

By f!4.39p and (I4.50H . setting + \l/ and taking into account Remark 14.31 we 

obtain exactly the weak formulation of problem fl4.llp - fj4.16p . □ 


Remark 4.5. Note that (14.281) is uniform with respect to 6. Moreover, we can obtain 
also estimates for WuJ t and S/yUfff uniformly in 5. Indeed, differentiating formally 
with respect to t problem (j4.18p ~ fj4.23p . multiplying equation (I4.18p (differentiated 
with respect to t) by ((uf t — and finally integrating by parts, we obtain, 

exploiting also the periodicity in time, 

i 



I Vuft + VyU S 't 


1,#|2 


da: dy dt < 7 


(4.51) 


0 12 Y 

where we used assumptions (12.71) . (12.121) and inequality (14.281) . Now, proceeding as 
in the proof of (I4.29p and (14.301) , we obtain 

1 



0 12 


Vufj 2 dxdt < y, 


(4.52) 


IIJ |V„uJ ;?| 2 dxdydt<y. 

0 12 Y 


(4.53) 


Therefore, passing to the limit for 5 —» 0 + , in (14.281) . (I4.52p and (14.53p . we obtain 
that the same estimates hold for 

This implies that (m # ,m 1,# ) belongs to C^([0,1]; H l {Q)) x C^([0,1]; L 2 (i7; X#(Y)). 

□ 


It remains to prove that any solution ( u,u l ) of the homogenized problem converges 
to (w#, w 1 ’#) as t —> 00 . This is the purpose of the next theorem. 

Theorem 4.6. Let{u,u l ) G L 2 (0, T); xL 2 (i?x (0, T); be the solution 

of problem (14.21) - (14.71) . Then, for t —y + 00 , (tqw 1 ) —> (u # ,u 1,# ) in the following 
sense: 


lim \\u(-,t) - u*(-,t)\\ H i {n) = 0 ; 

t —>+00 v ' 


(4.54) 


t —(“OO 


lim^ [||?x 1 (-, •, t) -u 1 ’ # (-,-,f)|| l 2 (dxV) + ||V 1 / it 1 (*, *, t) - V y u 1 ,# (-,-,t)||L 2 ( rt X y)] = 0 ; 

(4.55) 

(4.56) 


lim ||[u 1 ](-,-,t) - [« 1 ,# ](v,t)lU^Dxr) = 0 . 


































ASYMPTOTIC ... 


21 


Proof. Firstly we recall that, by m Lemma 5.2] which holds even in the present case, 
(u^u 1 ) G C°((0,T]-H\n)) x C°((0,T];L 2 (f2; Aj(y))) and [u 1 ] G C°((0, T]; L 2 (i? x 


ny 


As usual, let (r, r 1 ) := (w # — u,u 1,# — u 1 ), so that the pair (r, r 1 ) satisfies: 


cr (Vr + V^r 1 ) (V0 + V y <!>) dx d y d t 

o n y 

t t 

^ -pj^-^-[r 1 ][<F] dxdadt+a J J J [rl][$] dxdadt — 0 , VfG( 0 ,T), 

o n r o nr 

(4.57) 

where r = 0 on dfl x [0, T] in the trace sense, r 1 is periodic in Y and has zero 
mean value in Y for almost every (x,t) G i? x (0 ,T). Here 0 is any regular function 
depending on (x, f), with compact support in i? and $ is a any function depending 
on ( x,y,t ) which jumps across r, is zero when t = T and is regular elsewhere. 
Differentiating (14.5711 with respect to t, we get 





cr 


(Vr + Vyr 1 ) (V0 + V y $) dx d y + 



Q Y 


n r 


— fWl) [ r i] dx da 

[«!■#] -[U 1 ] L U J 


+ a [r]] [<f>] dx da = 0 . (4.58) 



n r 


Replacing (0,4*) with (r, r 1 ) in (14.5811 . we get 



ojVr + V yr 1 ] 2 dx d y + 



n y 


n r 


/(K’ # ])-/(K Dr 112, J 
-p-—- —- -r dxdcr 

[u 1 ’#] - [w 1 ] 


+ a J ][r]][r 1 } dxda = 0 . (4.59) 
n r 


As in Section^ equation (14.59[i implies that the function t H y a J n J r [r x (x, t )] 2 da dx 
is a positive, decreasing function of t, hence it tends to a limit value r 1 > 0 as t —> +oo. 
The value r 1 must be zero otherwise a J n J^r 1 ] 2 dadx > r 1 >0 for every t > 0. On 

the other hand, for t > 0 and setting T0 (t) |(x, y) G £2 x r : [r 1 ] 2 (x, y, t) < 2 a |r| |^| }; 

reasoning as in the proof of Theorem 13.21 it follows that 


a 


/ 


[r\x,y,t)] 2 dadx > y, 




Vt > 0. 


However, on T \ /0i, g(x,y,t ) := ^ > y > 0, where y is a suitable 

positive constant depending only on r 1 ,^, |T|, |i?| (this last result follows from the 
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assumptions (I2.8p ~ fj2.10p ). Hence, using (14.591) . we get 


cl I a 
dt 2 



[r 1 (x,i/, t)] 2 dadx < 


9 {x,y 1 t)[r l {x 1 y 1 t)} 2 da dx 


q r 


« r\r T i{t) 


< -x 


[r 1 (x, y, t)] 2 da dx < —— x < 0. (4.60) 

2a 


« r\r T i(t) 


Inequality (14.601) clearly contradicts the asymptotic convergence for t —> +oo of 
a Ur [r 1 ] 2 (a;, y, t) dadx to a positive number, hence 



lim a / / \r x {x, y, t)] 2 da dx = 0 , 

t— 4+oo ' ' 


(4.61) 


q r 


which is exactly (14.561) . Integrating (14.591) in [t, oo) and taking into account (I4.61D . 
we get 

+00 



cr|Vr + Vyr 1 ! 2 dxdydt < 


a 



[r{x, y, t)] 2 da da;, 


t a Y 

which implies 


n r 


+oo 


lim 

t —^-(-OO 



(4.62) 


(4.63) 


cr|Vr + S7 y r\ dx dy dt — 0 . 

t Y 17 

This last condition guarantees that for every positive rj there exists a t(rf) > 0, such 
that 

+oo 

r r 

cr|Vr + V,/| 2 dxdydt < rj , 

t n y 

which in turn implies that, for every n G AT", there exists a t n G (t + n, t + (n + 1)), 
such that 




a\Wr(x, t n ) + \7 y r l {x,y,t n )\ 2 dxdy < rj. 


(4.64) 


17 Y 


Hence, replacing (0,$) with (r t , r^) in (14.58p . we get 


a(X7r+ \7 y r 1 )('Vr t + X7 y rl)dxdy + J J g(x,y,t)[r 1 ][rl]dadx 

i7 r 



17 Y 


and 



cr(Vr + V y r 1 )(Vrj + S7 y r]) dxdy < 


+ aj = (4.65) 

i7 r 


———— [r 1 (a;, y, f)] 2 dcr da;. (4.66) 
2a 



17 Y 


i7 r 
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Moreover, integrating (14.661) in [t n ,t*], with t* G [t n ,t n + 2], we have 


sup 

t(E.[tn ,tn+2] 



Vr(x,t) + V y r 1 (x,y,t)\ 2 dxdy 


Y 


+ SU P I" [r l (x,y,t)] 2 d(rdx | , 

2 2a te[t„ )+00 ) 



Mn G N ; (4.67) 


12 r 


i.e., 


sup 

tG[t+l,+oo) 



cr 


|Vr(x,t) + V^r^x, y, t)\ 2 dxdy 


KSl Y 


“2 + a2 SUP 1“ [r x (a;, ?/, t)] 2 dcr ) . (4.68) 


te[t,+oo) 



si r 


Because of (I4.6ip the integral in the right-hand side of (I4.68P can be made smaller 
than ^ j provided t is chosen sufficiently large in dependence of rj. This means 

that 


sup 



a 


| Vr(x,f) + V y r 1 (x, y,t)\ 2 dxdy < rj. 


(4.69) 


Y 


Inequality (14.691) implies 


lim 

t —I - OO 



cr|Vr(x, t ) + V y r 1 (x, y , t)| 2 dxdy = 0 . 


(4.70) 


n y 


Now, working as done in (14 .31 p . we get 


lim / | Vr(x, t) | 2 dx dy = 0 ; and lim 

i t —^-)-oo J t —^-|-oo 

q by 



|V y r 1 (x,y,t)| 2 dxdy = 0. 


Finally, the previous results together with (14.611) and Poincare’s inequalities yield 

|r x (x, y, t)\ 2 dx = 0 , 



lim / |r(x, f)| 2 dx = 0 ; and lim 

t —^-j-oo J t-H1-00 _ 

n by 

which give (14.541) and (14.551) and conclude the proof. 


□ 


Remark 4.7. More in general, the previous procedure allows us to prove that solutions 
of (I4.2I) - (I4.7D having different initial data satisfying the assumptions stated at the 
beginning of this section but with the same boundary condition tend asymptotically 
one to the other (such convergence being exponential if / is coercive in the sense of 

(lull). □ 
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Remark 4.8. Observe that, thanks to previous remark, Theorem 13.21 implies unique¬ 
ness of the Deriodic solution (vft. of problem (14. 1 II) (14. 161) in^([0,l];^ 1 («))x 
C°([0,l];L 2 (f2;^(y))). □ 
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